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ABSTRACT 
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emphasis on . The oaner concludes by presenting some of 
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RESULTS AND DISCUSSION 


Governing Equations and Models 

Solutions to the Reynolds-averaged form of the Navier-Stokes equations written in mass- 
averaged (or Favre-averaged) variables are presented in the examples to follow. The goal of 
turbulence modeling is to relate the turbulent stresses and heat-fluxes to known mean-flow 
quantities such as velocity and temperature [5]. Equations can be derived to describe each of 
the Reynolds stresses [6]. If the Reynolds stresses and heat fluxes are related algebraically to 
the mean-flow variables, the corresponding models are called algebraic stress models 
(AS Ms). Eddy- viscosity models arc an important subclass of these models; they relate the 
Reynolds stresses to the strain rates in a manner identical with molecular stresses. This 
approach is generally referred to as the Boussinesq approximation. The models may be 
expressed in terms of an eddy-viscosity function and a turbulent Prandtl number. Compress- 
ibility is included by introducing mean density when mass-averaging is used. Aside from the 
use of mass-averaging instead of time-averaging, there is very little difference in form from 
the incompressible models. Eddy-viscosity noodels have been of primary interest for develop- 
ing the Reynolds-averaged Navier-Stokes codes to be discussed. Such interest results from 
considerations of ccHnputational efficiency and the limited improvement afforded by applying 
Reynolds stress models to complex flows. 

Compressibility and Averaging 

Many external aerodynamic flow problems involve speeds such that compressibility must 
be taken into account. Investigations by the author [7] and his colleagues, restricted to 
attached transonic and moderately supersonic flows, showed that compressibility effects 
could be accounted for by including the mean density in the mass-weighted model equations 
developed for incompressible flows [8]. This is illustrated in Fig. 1 where results for skin 
friction predicted by various models arc shown for a flat surface with and without pressure 
gradient. The predicted skin-friction variations without pressure gradient arc compared with 
the Van Driest prediction, which represents all available data to within 10%. Aside from 
showing that all nuxlels predict the trends of compressibility and pressure gradient with rea- 
sonable accuracy, it is concluded that the choice of mass- or time-averaging has no signifi- 
cant influence. The latter finding led to the decision to use mass-weighted variables in com- 
pressible Reynolds-averaged Navier-Stokes codes, since algorithms developed for laminar 
flows could be used and since modeling a large number of higher-oidcr correlations could be 
avoided. 

Transonic Flow Modeling 

Transonic flow computations are sensitive to the choice of a turbulence model. The high- 
speed flow leads to large displacement thicknesses, and shock waves can form that lead to 
separation. Each or both of these effects can influence the development of the entire flow 
field because of its elliptic character. Initial attempts to solve airfoil problems with Reynolds- 
averaged Navier-Stokes codes were disappointing (9J. In retrospect, the poor predictions 
resulted because the grids used were too coarse and because the turbulence models were 
inadequate. The situation has b^n ameliorated in both aspects. 

Fig. 2 shows predicted results using four eddy-viscosity models compared with the 
Bachalo- Johnson experimental data (case 861 1[9]). These recent computations were made 
with the same computer code using grids sufficiently fine to minimize numerical errors [10]. 
The Johnson-King [11] and Viegas-Rubesin [12] eddy-viscosity models do an excellent job 
of predicting this flow, which separates due to a strong shock wave. 

The deficiency of the Jones-Launder model is a result of the low-Reynolds-numbcr 
modeling terms and their integration near the surface. Viegas and Rubesin eliminated this 
deficiency by developing a wall-function model that includes the pressure-gradient effect and 
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Predicted results from the 0-Eq. models are not adequate, whereas the 2-Eq. results provide 
an improvement. In this model advance [5], the constant used to multiply the dilation term in 
the dissipation rate equation was modified to insure that the product of density and turbulent 
length scale remained constant during compression. Also, an upper limit was placed on the 
length scale in the eddy viscosity such that it never exceeded the Prandtl length scale in the 
near-wall region. This limit was necessary to control the maximum heating level at and near 
reattachment. 

Other studies of the influence of compressibility on models have concentrated on the 
problem of shear-layer development at supersonic speeds where significant compressibility 
effects arc known to be present. Vandromme [17] used Rubesin’s compressibility modifica- 
tion and modified the Jones-Launder model equations and in Fig, 8 his computations arc 
compared with shear-layer data (case 8500[9]). The spreading-rate from an unmodified 
Jones-Launder model does not show the decrease observed expcrimciitally. However, adding 
more compressibility terms to the model equations did give predictions that agreed with the 
data. More recently, Zeman has modified the dissipation rate equation to account for 
compressibility effects arising from shocklcts that form in supersonic shear layers, and he 
also predicts the proper spreading-rate [ 18]. The applicability of both of these model changes 
to the shock interaction problem discussed previously is being evaluated. 

CONCLUDING REMARKS 

Emphasis on model improvement for compressible aerodynamic flows using Reynolds- 
averag^ Navicr-Stokes c<^es has been directed toward eddy- viscosity formulations. As 
mentioned earlier, this stems from the requirement for computational efficiency and the per- 
ception that Reynolds stress models have not yet provided more accurate predictions. Com- 
putational efficiency will remain a requirement and, therefore, further improvement in eddy- 
viscosity modeling is needed. Some flows cannot be modeled with this approach, however; 
as a result, more accurate and computationally efficient Reynolds stress models are also 
needed. The challenge for the future is in the means used to achieve these needs. As the fol- 
lowing examples illustrate, I suggest that further work on the full- and large-eddy numerical 
simulations of turbulence will be the most fruitful. 

Most 2-Eq. model formulations are unable to model the near-wall behavior of turbulence, 
as some of our previous examples showed. Fig. 9 shows computations using various model 
formulations in the near-wall region of a turbulent channel flow compared with results from a 
full numerical simulation. Shih [19] (present model) was able to provide better formulations, 
using the simulation as guidance. Overall, his model predicts the proper behavior of the tur- 
bulent kinetic encr^, of other turbulence quantities, and of the eddy viscosity. This formula- 
tion has not been tried on any aerodynamic flow problems, however. 

Many aerodynamic flow computations will require Reynolds stress modeling to provide 
more accurate predictions. For three-dimensional flows, the cross-flow turbulent stress lags 
in response to changes in the cross-flow velocity gradient, and eddy-viscosity formulations 
may fail, depending on the severity and extent of the three dimensionality. Fig. 10 shows data 
[20] from a three-dimensional boundary layer experiment that illustrates the limitations of 
eddy- viscosity modeling. The cross flow, generated with a forward-spinning section of a 
cylinder aligned with the free stream, relaxes toward two-dimensional flow over the station- 
ary portion of the cylinder. The angles of the flow direction (0® coincides with the free 
stream), strain rate, and turbulent stress arc shown for a downstream location on the station- 
ary surface. Significant differences in angle between the stress and strain rate occur. The 
Rotta (1979) T-model (Fig. 10) predicts this difference between the stress and strain-rate 
directions, provided a judicious choice of reference frames is used (e.g., the spinning refer- 
ence frame). However, such a choice in this case depended on knowing the desired result and 
was not predictive. The Rodi (1976) algebraic stress model fails to predict the stress direction 
any better than an isotropic eddy-viscosity model, as indicated by the stress direction being 
nearly identical to the strain-rate direction. The Launder-Reece-Rodi (1975) model produces 
a better prediction, but is still not adequate. Numerical simulations for a flow analogous to 
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FIG. 1 Turbulent skin friction (from Ref. 7). (a) No pressure gradient, Rl - 10^; 
(b) pressure gradient, M«, = 2.3. 
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FIG. 2 Predicted pressure distributions over an axi symmetric circular arc bump compared 
with data. 
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FIG. 5 Pressure and skin fnction distributions on a two-dimensional compression comer 
deflected 20° (from Ref. 5). 
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FIG 7 Pressure and heat-transfer distribution on a two-dimensional compression comer 
deflected 15° at M« = 9.2 (from Ref. 5). 



FIG. 8 Spreading-rate versus Mach number (from Ref. 17). 
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FIG. 9 Near-wall turbulent kinetic energy and eddy-viscosity formulations for a channel 
compared with direct Navier-Stokcs simulations (DNS) (from Ref. 19). 



FIG. 10 Principal stress, strain rate, and flow direction in a three-dimensional boundary 
layer (from Ref. 20). 
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